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SUMMAEI 


The application of the Laplace transformation to the solution of 
^ the lateral and longitudinal stahility equations is presented. The 

espreaslons for the time history of the motion in response to a 
sinusoidal control motion are derived for the general case in which 
the Initial conditions, initial displacements and initial velocities, 
are assumed different from zero. Same illustrative examples of the 
application of the Laplace transform to ordinary linear differential 
eq.uations with constant coefficients and a numerical example of a 
specific prohlem are presented in appendixes. 


INTRODUCTION 


Recent developments in piloted and pilotless aircraft, eq.ulpped 
with automatic devices, have directed the attention of engineers to 
the theoretical investigation of dynamic longitudinal and lateral 
stability problems of aircraft designed for highr-speed and high- 
altitude flight. In the past, the dynamic stability Investigations 
were usually limited to the determination of Routh’s condition for 
stability and for the ceilculation of the roots of the character is tic 
stability equation to determine the damping of the modes of motion 
and the period of the oscillation. A more complete analysis of the 
problem requires the calculation of a time history of the airplane 
motion in response to a gust disturbance or in response to the 
application of the control surfaces. As the n^thods of classical 
analysis (references 1 and 2) proved to be inadequate for this 
purpose, new methods of operational mathematics, representing a 
more powerfril tool, were used. These methods are known today as the 
Heaviside operational calculus and the Laplace transformation. The 
application of the Heaviside operational calciilus to the calculation 
of airplane motions is discussed in references 3^ and 5 . However, 
the Laplace transformation is considered a more powerful method than 
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the Heaviside operational calculus hecause the Initial conditions of 
the problem, initial displacements and initial- velocities, are inher- 
ently taken into account by the Laplace transformation, whereas in 
the Heaviside operational cEilcTiliis, all initial conditions are zero. 

In this paper, the Laplace transformation is applied to both the 
longitudinal and lateral stability equations for the general case where 
the initial displacements and initial velocities were assumed different 
from zero. The operational equations obtained for this general case 
were then solved and the time hlsto 2 ry of the motion was obtained by the 
Heaviside expansion theorem and by the inversion theorem for_ Laplace 
transformation. The Laplace transformation is simple and effective. 

Its principles are easily understood and its technique quickly learned. 
It represents a further development in operational mathematics because 
it is a more powerfTil mathematical tool and because the difficulties 
and obscurities of- the work of Heaviside are avoided. 

A short historical sketch tracing the developnent of operational 
mathematics and its application to airplane dynamics is presented in 
appendix A. 

The author is Indebted to Mr, Leonard Stemf ield of-the Langley 
Stability Research. .Division, HACA, for information and collaboration 
he has contributed in connection with this paper. 


SYMBOLS 


c 

b 

S 

W 

m 

P 

Y 

t 


Sc 


chord, f eet ■ 
span, feet 

wing area, square feet 
weight, poimds 
mass, slugs 

density, slugs per cubic foot 
airspeed, feet per second 
time, seconds 

nondimens ional time parameter based on chord ^t ^ 
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St 

noncLlmerislonal time parameter based oa span ^ 


relatlTe densitj coefficient based on chord { 

\PSc 

) ' 


relative densitj coefficient based on span 


Dc 

differential operator with respect to Sq f j 

V<iSc/ 

1 


differential operator with respect to s-j, 


P 

operator in Laplace transformation 


\ 

root of stability equation 



radius of gyration about principal longitudinal axis, 
feet 


radius of gyration about principal lateral axis. 

feet 


radius of gyration about principal vertical axis 
feet 

9 


nondlmensional radius of gyration about principal 
longitudinal axis 


nondlmensional radius of gyration about principal 
lateral axis ^ky^/c^ 

^Zo 

nondlmensional radius of gyration about principal 
vertical axis ^kz^/b^ 

1] 

angle between principal longitudinal axis of inertia 
and flight path (fig. l), degrees 

6 

angle between reference axis and principal longl 
tudinal axis (fig. 1), degrees 

— 

a 

angle of attack (fig. l), degrees 
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flight— path angle 'betveerL path and. horizontal 
(fig. l), degrees 

attitude angle "between reference line and horizontal 
line, degrees (a + 7) 

deflection angles of- aileron, elevator, and rudder, 
degrees 

pitching angular velocity, radians per second (fl) 
angle of sideslip, radians 
azimuth angle, radians 

yawing angular velocity, radians per second (+) 
angle of "bank, radians 

rolling angular velocity, radians per second 
increment of forward velocity, feet per second 

nondlmenelonal increment“Of forward velocity 

nondimenslonal radius of gyration about longitudinal 
stahllity axis 

nondimenslonal radius of- gyration a"bout vertical 
sta"bility axis + Kj^^sin^r]^ 

nondimenslonal product— of— inertia parameter 


((%o^ - n cos ,1^ 


rolling-velocity parameter (helix angle generated 
"by wing tip in roll), radians 


pitching-velocity parameter, radians 
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1 

21, 7, Z 

Z 

Y 

Z 

L 

M 

N 

Cd = Cx 

Cx 

Cy 

^Z 

Cm. 

cz 

Cn 


yawing— Telocitj parameter, radians 


dynamic pressure 

.rectangular coordinates (fig. l) 
longitudinal force, pounds (fig. l) 
lateral force, pounds (fig. l) 
nonnal force, pounds (fig. l) 
rolling moment, foot-pounds (fig. l) 
pitching moment, foot-pounds (fig. l) 

yawing moment, foot-pounds (fig. l) 

* 

drag coe'fflclent 
lift coefficient 



longitudinal— force coefficient 


lateral— force coefficient 


normal— force coefficient 


pltchlng-moment coefficient 


rolling-moment coefficient 


yawing-moment coefficient 
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% = 




Scj 

he 


he 


'L. 


\27j 


“=^0 = 551 


hC^, 


"'Za = 
°ZDa 

^ze = 


% = 


ha 

hOr, 


hCz 

he 

\27 




^Z 

35e 


f - 


"“a 




BCm 

ha 


^Cjn 
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Cmg = 


^(i) 


W) 


K§) 


<§) 
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= 


he I 

^ Sa + 


Cn§5 


5 

ss;:®=‘ 


80^ 

he 


_l_ AA 


^ 5r 


Gl, G2, &3, Hi, Bq, H3 


^3 ^3 ^3 


functions of P on right side of oper- 
ational equations 


deteiminants 




S, C, P, s 


coefficients in fourth- degree characteristic 
equations 



abbreviated coefficients in operational 
equation 


Sn (H) • • • 833(H) 

^11 (h) ... h^^Cn) 


abbreviated functions of P in operational 
equation 


E 


residue 


The subscript o is used to indicate initial conditions, a bar 
is used to denote variables in the operational equations, and a dot is 
used to denote differentiation vith respect to time. 


AHALYSIS 


The purpose of this paper is to show how the longitudinal and 
lateral stability equations can be solved by the Laplace transformation. 
Thus no attempt is made to present a detailed discussion on the theory 
of Laplace transform, which can be found in references 6 and 7 and in the 
bibliography presented in appendix I of reference 6, but rather to 
present sufficient background of the theory to permit a clear under- 
standing of its application to this particular problem. 

If a function x(t), defined for all positive values of the 

variable t, is multiplied by e~^ and Integrated with respect- to t 
from zero to infinity, a new function x(P) of- the variable P is 
obtained; that is 
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x(P) = e ^x(t)cLt 

This operation on a function x(t) is called, the Laplace transformation 
of x(t) . The necessary eind sufficient conditions for the existence of 
the Laplace transform of a function x(t) are discussed in reference 6. 
Let 


^ ^ + . . . &nr-l ^ + anx = x(t) (l) 

dt°- dt^-^ it 

represen,x an ordinary linear differential eq.uation with constant 
coeff icl/ents b.-^, 8.23 a^, . . • a^j^. If D is substituted 

for — , for ^^3 and so forth, eq.uatlon (l) can he written in 

it d.t^ 

operational form 

(p^ + aj_D^^ + . . . + a^^x = x(t) (la) 

When t = 0, the following initial conditions ^e assumed: 


X = Xq 




dx 

dt 


^-1 


dt 


The Laplace . transformation of eq.uation (la), 
letter P for the operator, is 


with the use of the 
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(p^ + apl^"^ + . . . an_iP + an)3E = x(P) 

Term . 

corresponding 
to 

t 

+ (p^-^xo + P'^i + . . . Pxn_2 + Xj^i) 

tf^x 


+ ai(p^-2xo + + . . . Pxq^ 3 + x^g) 

D^-^x 


+ an-3(^^ + ?^1 + ^) 

d3x 


+ an_2(PXo + ^l) 

D^x 

• 

+ ®n— 1^ 

Dx (2) 

t 


The traneform i(P) for x(t) is taken frcan table I which presents 
some simple Laplace transforms. A more complete table of Laplace 
transforms is given in appendix HI of reference 6 and in appendix A 
of reference T. Appendix B shows two illustrative examples of the 
application of Laplace transform ttr ordi n a r y linear differential 
equations with constant coefficients. 


Longitudinal Motion 


The nondimenslonal linearized stability equations for longi- 
tudinal motion are given by NACA in the form: 

2u^Dc(a - 0) = 02^^53 + C^^^u* + | 


>( 3 ) 




'“9 2 “q 
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V The Laplace transformation is demonstrated for the case in which the 

elevator motion can bo simulated by the sine function 


6 q = Sitt sin asc (It) 

where 6 j^ is the amplitude. (In most cases 5^ ia assumed to be 1 .) 
Eearranglng and substituting equation (ll-) into equation ( 3 ) give: 

- %,)u' - (=Jc + I - (ozg + Iv„)8 

= 5 m® ®®C 

-CZu.»’ + [(2110 - I 0 zj 3 „)De - Oz^c. - + i 0 z^)d„ + Czg^ 

= Czo asc 


In order to illustrate the use of the Laplace transformation for 
a very general case, the only Initial condition assumed to be zero 
is 6 q^ = 0 ; that is, the deflection is measured from its trim 

position before the maneuver begins. For nil other parameters the 
initial conditions are assumed to be different from zero; thus 
the values are u^*, a^, 9^, and q^ at Sp = 0. The equations can 

then be written in general form, in which the four initial disturbances 
are combined with elevator motion, nil a specific problem some of the 
Initial conditions would probably be zero. For practical engineering 
piirposes, in fact, the most interesting cases are 

( 1 ) Disturbance only in angle of attack (due, for example, to a 
gust); elevator fixed; all other disturbances zero (Uq* = 0 q = ^ 


,1 'v 
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(2) Change In thrust, thus Uq* i- Oj elevator fixed j other 
disturhanoes zero = 0^ = q^Q = 0) 

(3) Dlsturhance caused hy elevator motion; other disturhanoes 

zero (uq* = = 0q = Io = 0) 

Each of these assumptions greatly simplifies the equations and 
shortens the computations, because many terms in equations developed 
for a general case vill vanish. 

The Laplace transformation of equation (3&) writtan. as 

follows : 


- Ox^,) u- - (Cj„ + I - (Oxg + I 0ij)8 


-Cz^,u’ + 


(^0 - I - 0z„]<^ - [(aio + I Cz^)p + 0^8 

! a 2 * (^° " 2 '^Zito)“o - (9*0 + I %) 8 o 

-(^,.'5’ - + 5 (a*oSy^ - 1 

- h ^7^ - I ^ SPoft2(pe„ + 10 ) 


(5) 


i V 


- i c,n_ao 


/ 




r 
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Equation ( 5 ) can be expreesed in a shorter form as 


where 


aiiu’ 

agiu* 

a3iu« 


+ + &13® = Gl 

+ aggot + ag^S = G-g 

+ a^gtt + ^ 23 ® “ ^3 




^11 = 

^12 = “(Cja I ^^jJ) 

^13 = -(^2e + I %^) 

^21 = -^Z^x 

=>22 = (a'c - I o^)p - 
^3 “ -[(s^c I =Z^)f Cz«] 

^31 = "^U» 

^32 " “(^ + 




^33 = 2^jjEy%2 _ i Q^p _ 
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(atpUo' - I - g CZq^o)(?^ 


(5l.= ^ 


(p2 + a2) 


-_8iM 

+ a' 


p2 j. 


gl(P) 


(P + la)(P - la) 


(7a) 


G2 =• 


aCZ5„^+ 


^2tic — ^ ^ZDa)“o ~ (^2^^ +“5 CZq^0oJ (P^ + a^) 


(p2 + a2) 


gg(P) 


(P + la)(P - la) 


(7b) 


G 3 = 


aCm5^5m+ f + <k>) ~ 


2 ^2^)a“'o 


- I %9^](r^ + “®) 


(p^ + a^) 




(P + ia)(P - ia) 


(7c) 


Now the eystem ( 5 ) or (5a) represents three simultaneous algebraic 

equations which can he solved for u*, a, and 0 hy the method of_ 
determinants . Thus 



A 


(8a) 




A 


(8b) 


(8c) 
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vhere tlie determinant 


^11 

&12 

^13 

^21 

®22 

&23 

^31 

^32 

^33 


The expansion of the determinant A results in a guartic equation 
in P 


A = AP^ + BP3 + cp2 + UP + E 

which generaJJy has two pairs of con^ilex conjugate roots^ namely 


^1,2 = -a i lb 
= -c i id 

Thus 

A = (P + a — ib)(P + a + ib)(P + c - id) (P + c + id) 
The coefficients of the quart ic (equation (lO)) 


A = 


15 


( 9 ) 


( 10 ) 


(10a) 


are 
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B = - ^Zu’^^Da) 

C = [^^^c(Cz^,Cmq^ + ~ ~ 

- Cn^% - CmD^% “ ~ •*• 

- %^Xu«^ZDa %°Zu«^XDa " %^Zu'^mDa ~ °Zg^inu>^XDa 
+ ^ZDa%%i) + 2UcKy2^Cx^,C2^ - ^Zui%) 

® = [l(^“a^u«% “ %%i^Za - ^^Xgi^^ZDa + ^mpa^Xu'^ ^ Sj^Z^i^Xa 

+ %C2^,Cj^ - Cj^C2^,Cjq^ - CxQ^z^.Cmpo, “ ~ ^ZQ^Tay;,i^Zj^ 

+ ^ZDaV«% ■" %. VXq) + 2|ic(CineCj^, + 

+ - Cj^C^g - - C^,Cj^)J 


E = %i^Cn^Czg - Cj^Cz^^ + Cz^,^CHjgCx^ - Cu^Cxg^ -H Cm^.^Cz^Oxg 
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The other deteiminants are 




®12 

^13 


Ai = 

&2 

®22 

823 

* (11) 

! 


^3 

^32 

^33 



^11 

Gl 

^13 


A2 = 

^21 

^2 

^23 

(12) 


831 


833 



^11 

ai 2 

% 


A 3 = 

®21 

322 

^2 

(13) 


^31 

®32 

‘'3 



When expanded the determinants can he written 

^ 2 . ~ (^) ~ O’P^ Tp (^) + 12(^5 (Ha) 

^ = — G’lfpi (^) + ^ 2 ^ 22 ^^^ ~~ ^3^23^^^ (l 2 a) 


= G-Lf^iCP) — G-2^32(^) + G-2f23(P) 


(13a) 
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where the minor e fj^j(P) ar'e 



®22 ®23 


ai2 &2.3 


®12 ^13 

fll = 


fl2 = 


^■13 = 



^32 ^33 


^32 ^3 


^22 ^23 



^21 ®23 


^11 ^3 


ail ^13 

^21 = 


f 22 = 


f23 = 



^31 ^33 


^31 ^33 


^21 ®23 



®21 ^22 


®11- ®12 


^ ®12 

f 31 = 

< 

^32 = 


f33 = 



^31 ^32 


®31 ®32 


^21 ®22 


After expanding, there result 



~ ^Da(t 


Cm00 Czj^ - 2^c) - Cina,(| 


+ 2n, 






( 14 ) 
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IT 

The other determi nants are 








ai2 

&13 



Ai = 

G 2 

822 

823 

1 

( 11 ) 


^3 

832 

^33 




ail 

Gl 

&13 



A 2 = 

®21 

^2 

823 


( 12 ) 


^31 

G3 

^33 




ail 

ai2 

Gi 1 



A 3 = 

a21 

822 

1 

Gg 


( 13 ) 


^31 

&32 

&3 



When e 2 j)anded the determinants 

can he written 



= Q-]f-[-] (P) — &pfTp(P) 

+ G-3f23(P) 

(Ha) 

A 2 = (^) + 

Gpfpp(P) — G 3 f 23 (P) 

( 12 a) 


A 3 = GjLf3l(^) ■" ^3^33^^^ 


(13a) 
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where the minora fj[^j(P) are 



^22 ^3 


^12 ®13 


ai2 &2.3 

^*11 = 

^32 ^33 

^12 = 

®32 ^33 

^13 = 

&22 &23 



^21 ®23 


^11 ^13 


®H ^13 

^21 = 


^22 = 


f23 = 



^31 ^33 


® 3l .^33 


^21 ^23 



^21 

&22 


^11 

®12 



^12 

II 

c — 1 

m 

^31 

^32 

^32 = 

®31 

®32 

^33 = 

^21 

^22 


After expanding, there result 


fn(r) 


[%(t “ ^c) ~ “ ^Za,^^cEY^ 

+ [Cingj^l Czjj^ - 2Pc) - Cz^ + 2\i^ + ; 

+ (^Za^me - CzgCin^j| (l4=) 


% 
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fl2(P) - |(-CX,^I^C%^)P3 

+ 5(%S * "^DaS - 

+ f Ja% “ %'’i%)j' (15) 

fl3(l') = |[cx4| % + Xc) + <=j^(x„ - 1 0zj^)]p2 

[%(l % ■ i °^) 

+ lfxDa% - 

(“la.% - %°2b.) ' (1®> 

f2l(P) . |(-4E2^,u^Ky2)p2 

- _°%'(2 % 2^o) - I Oz^.C^P 

-(V%'-%.S)j (17) 
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f22(P) 


- [t‘<=(2Cx,,%^ + S)]r^ 


* [!(%%• - 


+ 



(18) 


f2s(P) 


(%Cz^, - C2^,Czg)|. 


( 19 ) 


f3l(P) = 


2 

(%. s . - S .%)1 


1 n 

2 


) 


(20) 


^32(P) = 

- - I + I ^m^i^XDa)^ 

- - %i^D3a| 


(21) 
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, f33<« = 

+ [%.(| - ^Zj‘0 - I %.=I 

+ (%u.'=Za - %.%)| 


( 22 ) 


The solution of equations (8a), (8h), and (8c), which will result 
in a tima history of u* , a, and respectiyely, as a function 
of Sc, can he obtained frcsn. the Heayiside expansion theorem idien there 
are simple poles (reference 6). Uhls expansion theorem is an efficient 
method of findtog the inyerse Laplace transform of the quotient of two 

poljmomials — If, for example, 

I'(p) 


A 

where f (p) an(l_ F(p) are polynamlals with no common factors the 
degree of f(p) is lower than that of F(p), then for the case of 
simple poles and distinct roots 


= f- £2sL 

where are the linear and distinct roots of F(p) set equal to 

zero. The Heayiside expansion theorem is modified as indicated in 
reference 6 if any of the roots of F(p) = 0 are repeated linear 
factors. It is inportant to note that the expression for the Heayiside 
expansion theorem glyen here is different from the expression given 
in reference 4 hecauae of the different transforms of functions that 
are used in the Heaviside operational calculus and Laplace transfor- 
mation. However, if a problem is consistently followed through by 
either one of these two operational methods, identical solutions will 
be obtained. 

The application of the inversion theorem of Laplace transformation 
to the solution of equations (8a), (8b), and (8c) by ccanputing residues 
is shewn in appendix C. 


f(p) 

bTpT 
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Lateral Motion 

The nondimensional linearized HACA standard eq.uations of motion 

are: 


Sideslip; 

+ V) = +i +Ci,i<+ (Cl tan r)t+ | Or^V 


EoU: 




>(23)- 


Yaw: 


2n^(l2%2^ + EjzLb20) =C^^6 + C^P + | + | 


It is Import^t to note that from the standpoint- of mechanics, 

Kj 2 , should he defined as Kxz = “(^Zq^ “ '*1 ®°® ''1 & right- 

hand system of axes. However, the definition of- Kxv. presented in 
the. symhol list is used in the paper to conform with recent^ NACA standard 
equations of motion. 

The Laplace transformation is demonstrated for the case in which 
the control-surface motion can he simulated hy the sine function 

5 = sin as -^3 (24) 

where is the amplitude, (in most cases &m is assumed to he 1.) 

Rearranging and suhstitutlng for 5 give; 
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~ (2 ■*’ “ 2 ~ ^ - ^Tg^ ®®b 

“® 2 pP + - I “ (2 ~ 2 M.-gS:j 2 ;Db^^* = sin asg , I ^g^a) 


~ (2 ~ 2^tEizDb®)f^ + \ <kpl>b^+ = CngSjj sin aa^ 

V 

If the initial coiditicjns (for Sg = O) are 0^, ijr^, p^, and r^ and the trim position 
is assumed as zero (thus '5 q ° O), the Laplace transfccmaticii (with P substituted for p 
to aToid confusion with angular Tolooitf) for equation ( 23 a) can he written as follows A 


^ 2 ngP - Oyp)p - Oj^P + Ol)0 + (2ni^ - i Cj^P - Cl tan r)+ = OjgB^ 
+ 2u^jPo - 2 °Tpi^o + (^IJ'h “ 2 °Tr)+o 

-ClpF + (2uhEl2p2 - 1 Cz^p)? - Cz^P - 2 m.z>K:x^)v = CZgSm ^ 2 “^ 

+ 2 vxgKj 2 ^pj(^ + Po) - I Olpi^o - I Czr+o + + ^o) 

-PnpP - (1 %P - 2w,E:x^)^ + (a^hEZ^P^ - 1 On^p)t = Qng5a g 

- ^ Pnp0o + 2jij,^XZ^P0o + Po) + + ^o) " 2 *^r*o 


Equations ( 25 ) can he ezpressed in shorter form as 


tn P + hip^ + hLg''!' = Ej 

bglP + hgg? + hggt = 

hj^P + ^220 + h 33 ^ = H3 


V ( 25 a) 


^or practical engineering purposes a singdlfied case is of interest, namely 
respanse to a horizontal gust + 0 , -rfiile Pq = 0o = Po = ^o = 6 = 0. 
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where 


- Cyp 

tl2 = CYpP + Cl) 

^13 = (2^^bP - I Cy^P - Cl tan r) 
h2i = -Cjp 

b22 = _ 1 cj^p^ 

>>23 = -(I 

I331 = -Cnp 

»32 = -(I - 2 i*b%zr^) 

»33 " - I “nrl) 


(26) 


% - 




_2tibPo - I Cyp0o - (I Cy^ - 2Hb)to](p2 + a2) 


(p2 + a2) 


hl(P-) 

P2 4. a2 

hl(P) 

(P + ia)(P - ia) 


(27a) 
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Clg5ma+ +Po) -| 2n^jZ^Piiro + r„)] (p2 + a2) 

p2 + a2 


h2(P) 

P2 + a2 

= ^2(P) (2Jh) 

(P + la) (P - la) 


M = 


Cns^a+[-| Cr^0o + 2UTfe(P0Q+Po) + 2tx^2(Pto + rp) - | (pS + a^) 


+ ar 


ti3(P) 
p2 + a2 

h3(P) 

(P + la) (P - la) 


(27c) 


Now equations (25) or equations (25a) represent tkree simultaneous 

algetralc eq.uatlons which will he solTed. for P, 0, and ■'Jr hy the 
method of determinants 

( 28 a) 

A 

0 = ^ (28b) 

A 

A 


(28c) 
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where the determinants are 


'll 

bi2 

bi3 


'21 

^22 

^23 

(29) 

’31 

^32 

^33 


Si 

bi2 

bi3 


Eq 

b22 

^23 

(30) 


^32 

^33 


^11 

Si 

bi3 


b2i 

^2 

^23 

(31) 

b^i 

% 

^33 


^11 

bi2 

%; 


^21 

b22 

B2 

(32) 

”31 

^32 




If the values of equations (26) are substituted into equation (29) and 
the determinant A is es^janded, a quartic equation is obtained 


A = AP^ + BP3 + CP2 + DP + E = 0 


( 33 ) 
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which, generally has a pair of conLplex conjugate roots and. two reeLL 
roots 


P ]_^2 = -a ±. ih (Batch.— roll osclllaticoi) 
P 3 = -c (spiral mode) 

= -d (rolling subsidence) 


Thus 

A = (P + a - lb)(P + a + ib) (P + c)(P + d) (33a) 

The coefficients of e(iuatlon (33) are 

A = 8Uh3(kj%^2 _ 

B = + On;^) + Iz^Jp - 

0 - - 0^0^) + - %0lp) 

- Ssz^OYpOj^ + OYpOnp + ‘‘MbOlp - 0Tj,Cjp - Ol^Onp) 

- - %og] 
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D = |-2^^CLjtan ~ ~ ^Z^npJ 

+ + "2p(%% - %S) 

+ Cnp^Cy^Cip - CypCj^^J + hb(^2p% " C^jpCip^^ 

E = ^jtan 7|^CnpC2p - CipCn^^ + CipCn^ - CnpCj^J 

Kie development of the determinants (30)j (3l), aJid (32) gives 

Ai = HiFii(P) - HgFxgCP) + H3 Fi 3(P) (30a) 

^2 = -"HtFqt (P) + ” 23^23(^5 (31a) 

A3 = HiF 3 i(P) - B^ 32(P) + H^33(P) (32a) 


where the minors are 



^22 

^23 

Fi 2 = 

hi 2 

hl 3 

F13 = 

hl 2 

^13 

Fii = 

h 32 

•V- 



^33 


^32 

^33 

^22 

^23 


h21 

h 23 


^11 

hl 3 

F23 = 

fell 

bl 3 

F21 = 

IJ3I 

^33 

F22 = 

^31 

133 

tgi 

^23 


h 2 i 

^22 


I’ll 

hi 2 


^11 

’*’12 

II 

1 — 1 

m 

^31 

132 

F32 = 

^31 

^32 

^33 " 

h 2 i 

^22 
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The es^aasloas yield 

Fii(p) = ■'■ 

- %Z(% + Cj^|p3 + - Ci^C^)p2 (31^) 

^12^^^ = “ ^Yr) “ 

+ - CY/np) - ^^(2Kz2CL - ^JzPl tan 7 - 

+ "^(pn-r “ 


+ CYj,) - ^h)^XZCyp]p^ 

+ -ft(C2rCYp - CypClp) + U^|2CL(Ex^tan 7 - Kxz) + 


^P2 


+ - CZptan 7)P 


(36) 


P2l(P) = -2V^b(^^Zp - E2zCnp)p2 + |(CzpCn^ - C^pOgP (37) 
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I'22(P) = + C^)p2 

+ [|(°nj.% - “Y/np) + 2l*b%]r' - =1 tan 7 C^X (38) 

Y23(r) = + l*b(-2H:izOYp - 

* |(°ll.°Yp “ '’Yr'^Ip ■’■ - Oj, tan 7 Cj^ (39) 

r3l(P) . [2Mi(-«zz0lp + El^np)?® + |(“npCjp - CnpCjj,)p] (40) 

P3a(P) - [w%z^'^ + Mt(-2K:3EOYp - C^)p2 + |(c„j,0yp - CnpOyj,)P 


P33(P) = _ i^(2K^%^ + Cjj,)p2 + |(CYpCj^ - %Cjp)P 



(42) 


As Indicated in the section entitled "Longitudinal Motion" the 
solution of equations (28a), (28h), and (28c), which will result in 
a time history of p, 0, and t, respectively, as a function of Sij, 

can be obtained ffam the Heaviside expansion theorem (reference 6) 
or by computing residues as shown in appendix 0. 



NACA TN 2002 


31 


General Eemarks 

The method, presented, fji this paper can also he applied, to cases in 
vhich the airplane is equipped, with automatic pilots acting on elevator, 
rudder, and ailerons. Each automatic pilot is characterized hy addi- 
tional equations of motion. Thus, generally speaking, there are four 
simultaneous equations for longitudinal motion and five simultaneous 
equations for lateral motion (two automatic pilots) . As before, the 
Laplace transformation is applied to these equations and the problem is 
treated according to the method indicated in this paper. 

Some transfoims of simulating functions for control— surface motion 
are presented in appendix D and simplified methods of computation of 
Laplace transformation are given in appendix E. 


CONCLUDUKJ BEMAEKB 


The application of the Laplace transformation to the solution 
of the lateral and longitudinal stability equations has been presented. 
The esqjressions for the time history of the motion in response to a 
sinusoidal control motion were derived for the general case in which 
the initial conditions, initial displacements and initial velocities, 
were assumed different from zero. 


Eyan Aeronautical Company 

Lindbergh Field, San Diego, Calif., July 22, 19^9 
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APPENDIX A 
KESTOEICAL SKETCH 


A short historical sketch on the development of the operational 
calculus and Its application to airplane dynamics is presented. 

The fundamentals for the theory of «rnan oscillations about a 
steady state of motion were developed- in I877 by South (references 8 
and 9). Then as early as 1903 Bryan applied the mathematical equations 
of motion of a rigid body to the disturbed moticn of an airplane 
(reference lO). In the following years the mathematical theory remained 
fundamentally in the form proposed by Bryan, but the method of appli- 
cation was changed as the reaiilt of- the development of experimental 
research by the NACA. 

During those years many scientists were working on the problems of 
dynamic stablll-ty, not only in the United States but also in Great 
Britain, France, Belgium, Germany, and other countries. In 1927 the 
eq_uatlons of motion were first expressed in dimensionless form 
by Glauert (reference 11 ). Jones, Bairstow, Zimmerman, and Millikan 
(references 1, 2, 12, 13, and l^t-) also dealt with-dynamic stability 
and their work is well-known to the average engineer in this comtry. 

The need for a meeins of describing the response of the system 
(mathematically similar to the system iised herein) to the applied 
disturbance was realized by electrical engineers many years ago. 

In 1899 Heaviside, Impelled by this need, contributed a significant- 
development. In his electromagnetic theory he originally devised 
his operational calculus for the solution of ordinary linear differ- 
ential eq.uations with constant coefficients and some of the partial 
differential equations of applied mathematics. The principles of 
this method are illiistrated in reference 15* 

The significance of Heaviside's contributions were not recognized 
in his lifetime because of the inadequacy of- the mathematical treatment 
and the obscurity of-his papers. Bromwich, making use of the theory of 
functions of a complex variable, expleiined and established the validity 
of Heaviside's methods. Bromwich's method consisted of finding the 
solution of a given differential equation, with initial and boimdary 
conditions, in the form of a complex Integral over a suitable path; 
the choice of the Integrand and contour is sometimes difficult. Further 
reseeirch by Carson, Carslaw and Jaeger, and Doetsch (references 16 
to 18) resulted in the application of the Laplace transformation to 
the differential equation. Finally, Doetsch recognized fully the 
value of the "inversion theorem" for the Laplace transformation. Thus 
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■fche Laplace ■fcransforma'tlon is an important step forward in operational 
mathematics. A complete treatment of the subject of Laplace trans- 
formation can be found in references 6, 7, and I 7 . For some time now 
it has been recognized that by applying the Laplace transform a better 
substitute for HeaTiside operational methods can be obtained. 

Among the early attempts to apply operational calculus to the 
problems of stability and control was a very fundamental work well- 
known in this country (reference i|-) . This paper deals with lateral 
motion and applies Heaviside operational calculus. Later several 
papers were written on the dynamic response of the aircraft which also 
made use of the Heaviside method. Some dealt with tail load variations 
due to elevator motion (for example, see reference 19) . Others dealt 
with stability with free controls (reference 20), stick forces in 
maneuvers (reference 2l) , and the behavior of the airplane equipped 
with automatic control (reference 22) . 
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APPENDIX B 

ILLDSTEATIVE EXAMPLES OF APPLICATION OF LAPLACE TEANBPORMATION 

Example I 


Example I illuBtrates the application of Laplace transform, to the 
equation ^ 


(d 2 - 3D + 2)x = eat (t < O) 


The initial conditions when t = 0 are 


X = Xo 
Dx = x^ 


Table I (transform 3) shows that the transform of 


x(t) =e&t 


is 


x(p) 


1 

(p - a) 


Thus the Laplace transformation of the given equation is (n = 2 when 
equation (2) is applied) 


(p2 - 3p + 2)x = + (pxo + xi) - 3 Xq 
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Example II 

Example II applies the Laplace transform to the equation 
(d3 - 2D2 + d)x = il- (t > 0) 

The Initial conditions vhen t = 0 are 

Xq = 1 

xi = 2 
X2 = — 2 

Table I (transform l) gives the transform for »• 

x(t) = k 
as 

x(p) = I 

If the rules of equation (2) are applied, the transform can he 
written (n = 3) as 

(p3 - 2p2 + p)5E = ^ + (p2xo + pxj^ + Xg) - 2 (px^ + x^) + x^ 

= I + (p^ + 2p - 2 ) - 2(p + 2) +1 
_ p3 — ^p + !<■ 


P 
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APPENDIX C 

THE mVEESIOH THEOREM OF LAPLACE TRANSFORMATION 


By means of-^this theorem x(t) can he obtained from its 
transform 3c(P) . If 


then 


x(P) 



0 — Ptx(t)dt' — 


R(P) > 0 


x(t) 


1 

2i« 


J p7+ioo 

x(X)e^^dX 

7-i« 


(Cl) 


where 7 is a constant greater than the real part' of all singularities 
of 3c(X) and 


~17+i« 

7-ioo 


lim 

W 00 



The path (7 — ioo, 7 + ioo) may he replaced hy a circle C containing 
all the poles of the integrand. Then x(t) is equal to 2ni times 
the sum of- residues at these poles. The method of evaluating the 
residues is shown at the end of this appendix. 


Example I 


J - P3 P - 4 
p2 - 2P-- 3 


Let - 
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The inversion theorem, with X auhstituted for P in the eq.uation, 
gives 


z = -3- ^ e^^dX 
2i« Jy_ico X2 + 2X - 3 

The denominator has two roots —3 and 1. The residues must he evaluated 
at two simple poles at X = 1 and at X = — 3 and then summed in order 
to obtain x. 


Example 2 

Consider the simultaneous equations 


(3p + 2)x + Py = i 


+ (4p + 3)y = 0 


which yield 


+ 3 

P(P + 1)(UP + 6) 


(UP + 6)(P + 1) 


Application of the inversion theorem of Laplace transformation gives the 
solution 


X 


1 

2lJt 


37+1“ 
7-1 « 


i^X + 3 

X(X + 1)(UX + 6) 


e^’^dX 
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y 


J7-I00 (llX + 6)(X + 1) 


The sums of residues on each, pole for 2; and y give the solution as 
a function of time. 

The inversion theorem for Laplace transformation is now applied 
to equations (8a), (8h), and (8c). If A. is suhstituted for P, new 
determinants A, A^^, Ag, and A^ are obtained from A, Aj^, S2, and 

A time history of u* , a, and d as a function of the parameter Bq 
is obtained when the inversion theorem is applied to these equations 
(the path angle 7 is then also deteMnined as 0 = a + 7) : 


u‘ 


1 

21it 


a 


1 

2iit 


p7+lco^ 

j7-ioo ^ 

U7-i«> ^ 




e^"°dX 


(C2) 


(C3) 


0 


1 

21s 


lr-i» ^ 


(04) 


In order to evaluate the integrals (C2), (C3), and (C4), the values of 
all residues for each integral must be found and summed. This procedure 
is demonstrated bn equation (C2) . The parameter \ is substituted 
for P, equation (lOa) is substituted for S, and equation (lla) is 
substituted for The auxiliary substitutions for G]_(X), G-2(X), 

and G^(\) are obtained from equations (ja)j (Tb), and (7c). The 
values for fll(X), fj^2(^)j fi3(^) are obtained from equa- 

tions (1^), (15) i and (16) and \ is again substituted for P. Then 
equation (C2) can be written 





% I 


u‘ = 


1 

2ijc 


">7+1“ 


ioo (X + la)(^X — la)(X + 1) — ic)(X + b + 1g)(X + d — le)(X + d + ie) 


(C5) 


The residues B must he evaluated for six simple poles; at pole at X =» — ia. 

Eg at pole at X = la. Eg at pole at X = -h + lo, at pole at X « -h — Ic, 

E^ at i>ole at X = -d + le, and Eg at pole at X = -d — le. Finally u* = ZE. With 
ttie use of the relations 


sin I = 



cos X 


elx + e-^ 
2 


the u‘ terms can he collected and expressed as sine and cosine terms, as illustrated in. the 
numerical example. If any function other than the sine function vere used to simulate elevator 
motion, there would he no changes in equations (lla), (l2a), and (l3a) Qfor and Ag^ 

hut the expressions for the functions Gr-^, Gg, G^ (equations (Ta)^ (7t), and (7c)) would he 
different. (See appendix D. ) 


LO 
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Numerical Example 


1 9^^c(x 2 + I.4773l2gx- ’i- i3^^A908846)d>. 

Jy— lob ~ ^l)\^ “■ ^ 2 )(^ “ ^3) ~ H) 


The roots of the quartlc are 

^1,2 = -0.739517 ± 11.56351 
^3,4 = -0.016582 t 0.044189561 

Suhstitutlng these values for X^t X^, and. gives 

1 e^°c (xg + l.X7731gg>- + l34.4908846)dX ______ 

2*1 Jpc-i* (>■ + 0.739517 - II.56351K1 + 0.739517 + U.56351)(1 + 0.016582 - O.OiAl8956l)(X + O.Ol65^ + 0.(Ml8956i) 
In this example there are four simple poles. First the resioLue Is ccmputed, at the pole 

Xi = -0.739517 + 11.56351 


Then 


% = 


q>-1Sc(j^^2 + 1.4773122X3^ + 134.4908846) 

(Xj + 0.739517 + u. 56351) (Xi + 0.016582 - 0.044189561) (Xi + 0.016582 + 0.044189561) 
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Sutetituting the value for In the last sq.uation yields 

^ 0^isc[_i33. 1676468 - 17.102809661 + 1.4 t 73122(-0. 739517 + 11.56351) + l3l^.l^.9088l|■6’ 

^ “ -3080.28381 + 386.65737 

- 0.2307403 - 0.019911) 

-3080.28381 + 386.65737 

With the use of the algehralc manipulation of complei numhers 


E « gXiSc (0.2307403 - 0.019911) (3080.28381 + 386.65737) 

^ (-3080.28381 + 386.65737)(3080.2838 i + 386.65737) 

= e>-lBc (0.1505458 + 0.703047331) 

9637.653 

= (O.OOOOI562 + 0.00007294l)e^~^’^^^^^^‘^^*^^^^^^®c 
At pole X2 = -0. 739517 — 11.56351 the sign of imaginary part 1 b changed, thus 

B2= (0.00001562 - 0.000072941)6^“®*^^^^^^“^*^^^^^®® 

At pole X3 = -0.016582 + 0.044189561 

^ gX3Bc (X33 + 1.4773122X3 + 134.4908846) 

3 ® (X3 + 0.739517 - 11.56351) (>.3 + 0.739517 + 11.56351) (I3 + 0.0165824 + 0.044189561) 
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Sul)stitirblng for yields 


Bo = 0‘ 


X3SC [-0-O0I6775 - O.OOlhe^^ + 1.47T312g(-Q.Ql6582 + 0.044l8956l) + 13^^.4908846] 


11.8173841 - 0.0056074823 


PoUowing the same procedure as that used for the X]_ root gives 

Bo = (134.46471 + O.063816281) (-11.8173841 - 0.0056074823) 

(11.8173841 - 0.0056074823) (-11.8173841 - 0.0056074823) 


X3SC (-^.( 


-139.650596 


= (0.0000009546 - 11.378551)0^”®*®^^^^'®^^®^^^^^®° 

At pole Xjj^ =» —0. 016582 — 0.044189561 the sign of the Imaginary part Is changed, thus 


E4 = (0.0000009546 + 11,378551)6^"®*®^^^“®*®^^®^^^^^®° 
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Now u’ = + E2 + Substituting tbs values for R^^, E2 j‘ 

and R]^ yields 

u = (0.00001562 + 0.000072948i)e^“^*'^^^5^’^‘^^*^351)sc 
+ (0.00001562 - 0. 0000729481) e^“°* '^39517-11 *56351) Sc 
+ (0.0000009546 - ii.3Y355i)e(-0-01^582+0.044l8956i)sc 
+ (0.0000009546 + 11 . 378551 ) 

u» = e“°*'^^^^^'^®°(o.00001562e^*^^^^^®° + 0.00001562e“^^’^^^^^®° 

+ 0.000072948le^^‘^^^^^®o - 0.000072948le"^*^^^^®°^ 

+ e"° • 016582s c (q ^ O 000009546e“^ * 044l8956isc 

+ 0.0000009546e"°‘°^^®5^^®° - u, 3785510® 

+ 11. 3785510“® •®^^®556iBc) 

Wltb the use of the relations 

el^ + e—l*- = 2 COB x 

e^ — e“^ = 21 sin x 
the value of u’ ceui be expressed as 

u‘ (0.00003124 cos H.5635sc-afOOOl45896Binll.5635ec) 

+ e“®‘®^^^®^®° (0.0000019092 ops 0.044 i 8956 bc + 22.75710 sin 0. 04418956s c) 

The roots must be computed very exactly to several decimals; other- 
wise the computation by the method,, of residues does not check to zero. 
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Evaluation of Besldues 

(a) Simple pole: 

The residue at a simple pole of the function f(z) is 

E = 11m (z — a)f(z) 
z >a 

(h) Multiple pole: 

Let 

g(z) = (z - a)“f(z) 

Then the residue at the pole of mth order is 


(C 6 ) 


(C?) 




(m- 1 )! z-»a 

g(z)l 


1 

^m-l 

(m - 1 ) 1 



z=a 


(C 8 ) 


Example 1 (Simple poles).— Let 


X 


1 

2 irt 


07+loo 

J 7-1 00 


X3 + X - 4 

(X 2 _ 2 X + 2 ) (x 2 + 2 X - 3) 


e^^dX 


The denominator has four roots. There are four simple poles (m = l) 
at X]_ = —3, X 2 = Ij and 13^4 = 1 ±. i. The computation of the residue 

at a simple pole X = — 3 is illustrated. If- the term under the integral 
sign is called f(X)dX, then according to eq,uation (C7) 


g(X) = (X - l)f(X) 
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Tims the residue is 


(\ - l)(x3 + X - 

U2 _ 2X + 2) (X - 1) (X + 3)J^^_3 


(-gT - 3 - 
(9 + 6 + 2)(-3 - 1) 


= -r=2t_ e-3<^ = i 6-3* 

-4 X 17 2 

Example 2 ( Double pole) Let 


(t) = -L. (xxq + XI + 4xq) 

2i« U7-lco (X + 2)^ 


e^''=dX 


There is now a double pole (m = 2) at X = a = According to 
equation (Cy) 

g(X) = (X + 2)2f(x) 

which, when substituted into equation (C8), yields 


R = i ^ + 2)^(>^3 Co + XI + 4Xq) 

(X + 2)2 


= + t(Xxo + + 4xo)e^'^ 


X=-2 


= Xoe”2t + t(-2xo + + 4xo)e" 

= [xq + (2Xq + xi)t]e“2^ 
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Example 3 (Triple pole) Let 


X = 


^ >; e^’tdx 

2in ,L , . , ^ «v3 


(x2 + a2)^ 


There are triple poles at X = ia and at X = — ia. If X 
then equation (C7) gives 


g(X) = (X - la)3 L 


At 


(x^ + a^)^ 

which, substituted into equation (C8), gives 




(X - ia)\e^'^ 


gx + ia)(X - ia)]- 


X=ia 


l6a^ 


* a)' 


dat 


Similarly, at X = — ia 


Ep - -V 

16 ?^ a/ 


iat 


X = E2_ + Eg 


With the use of the relations 


sin X = 


gix _ Q-ix 

2i 


cos X = 


e^ + e ^ 


I 


= ia. 
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APEEaroiz D 

TEAH5F0EMS OF SIMULATIMj FIINCTIOIJS FOE COHIEOL-SURFACE MOTION 


In order to include the control-HSurface deflection as a function 
of time ^or parameter Sq and s-^) in. the eq.uatlons of motion, the 

assumed motion must he simulated (for simplicity) hy some simple knxjwn 
functions. 

Some examples are given vhich could he satisfactory in many 
practical cases. The functions are given here for five examples and 
their transforms are given in table II. 

(l) Step function: 



(F = Oj t < 0) 
(F = Fq; t > 0) 


(la) Straight line: 



(F = mt + h) 



4 . 
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( 5 ) 



It is often convenient to get the result for unity of control- 
surface deflection (say one radian); then tlie result for any arhitrary 
deflection can he readily computed. 
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APEEINDIZ E 

SIMPLIFIED METBDDS OF COMPOTATIOIT 


The computation of the time history of any parameter can he 
shortened in some cases if the foim of the solution and also the 
value of the function and its first, second, and third derivatives 
at a given time are known. It is preferable to find these boundary 
conditions at time t = 0 (initial values). On the assumption that 
the Inversion theorem for Laplace transformation gives a parameter y 
in the form 


y = 



7+l«> 


f(X)e^‘‘^dX 


y_ico (X - Xi)(X - XgXX - X3)(X - 


(El) 


this method will be applied to the longitudinal and lateral motions. 

Longitudinal motion .— The denominator (of the stability eq.uation) 
has five roots with one root X5 = 0. The other four roots are 
freq.uently two conjugate complez pairs: 

Xj^ ^ ^ CL i. ip 

X3^i, = 5 ± ir 

If eq.uatlon (l) is assumed to give the solution in the form 

y = Cq + e°'^(Cx sin pt + C2 cos Pt) + 0^(03 sin 5t + C4 cos 5t) . . . 

(E2) 

then Cq = lim y can be evaluated. Differentiation of eq.uatlon (E2) 
t-^“ 

gives the derivatives y, y, and y*. 

If the initial values of these functions are known, that is, yg, 

7o> y*o^ y*oJ t = 0 the constants C2, C3, and Cij_ of 

. I - , . I , « ■■ ■ _ I II .1 — — — ' 

^The author is Indebted for this suggestion to Mr. J. M. Debevoise. 
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equation (l2) can be obtained by solving foirr simultaneous algebraic 
eq.uations which can be written (for t = O) as follows: 


= Co + Cg + C]|^ • • • (E3) 

yQ = + otCg + + 7Ci|_ ... (E4) 

y-Q = 2a3Ci + (a^ - p2)cg + 276C3 + { 7 ^ - 62)C]^ . . . (E?) 

y‘o = ^ 3 a 2 p - p 3 )c^ + (o 3 _ 3 ap 2 )c 2 + (3725 - 53)03 

+ (y3 _ 3752)02 ^ . . . (e 6) 


The value of y* can also be obtained by plotting y against t. 


y = e 


|j[a2 — p2)cj^ _ 2apC^sin 3t + j^2apC;i^ + (a^ — p2)Q^ qqq 
^ " ^^^^3 “■ 278Ci^sin 6 t + [ 278 C 3 + (/^ - 62 )Ci^cos 6 t^ 


Lateral motion .— For the case of lateral motion the equation has 
one root = 0, one conjugate complex pair, and two real roots 


^1,2 = a ± 

X3 = 7 
X4 = 6 

The time history of any parameter y can be written 

y = Cq + e^^(C2^ sin Pt + C2 cos Pt) + s^C3 + (E7) 
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and., as "before, 0© = 13m y. The initial values of 7 q, y^, yn* 
an d yQ for t = 0 are as stoned to "be known. Differentiating 

eq,uation (EJ) and sutstituting for t = 0 give four simultaneous 
algetraic equations from which, the constants C 2 , and C)j. 

of equation (E2) can he determined: 

Jq = Cq + C 2 + + C]j^ ... (e8) 

yQ = PC]^ + C 1 C 2 + + ^1). • • • (®9) 

y^j = 2opCi + (a2 - 3^)02 + 7 ^^ + 

■fo = (3a2p - p3)c^ + (a3 _ 30 ^ 2 ) 0 ^ + 73c^ + . . . (eh) 

If desired the value of *y‘ can he obtained hy plotting y against 
time 


y = 6°* |[(a2 _ p2)c^ _ 2opc^sin pt + [2aPCi + (a^ - p2)cj cos pt 

+ + e®''^52ci^ 


In the cases in which there are simple poles, the .method presented 
in reference 7 can he used. This irethod is Illustrated briefly. It 
has been seen that any parameter y of the equations of motion (such 
as a, 0, ilf, and so forth) es^ressed as a function of the_ 
operator P can he written as a ratio of two determinants A(P) 
and B(P) which are polynomials in P 


y(P) 


a(p) 

B(P) 


Thus the characteristic equation is 


1(P) = 0 
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This eq.iiatlon is a polyncaaial in P, the highest pover of P "being q_; 
thus it has q. roots: X^_, Xg, . . . X^. The inversion theorem gives 


y(t) 


A(t) ^ ^ 
B(t) 


(t >0) 


wliere 


B* = ^ B(P) 


For multiple poles and special cases, see reference 7 (pp. 152 

to 169 ) . 


« 
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TABLE I 


SIMPLE LAELAOE TRABSFOBMS 
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TABLE n 

TEAHSFCERMS OF SIMOLATIHG FIMCTIOIB FOR COWTROL-^STIRPACE MOTION 


Type 

Control-aurface 
motion determining 
x(t) 

Transform for 
initial condition 
5o = 0 

x(P) = e“^r;(t)dt 

* 

Transform for 
assumption 
Fo = 1 
x(P) 

1 

FQ(t) step 
function 

^*0 

P 

1 

P 

la 

mt + t 

(m + TdP) /p2 

For t = 0, m/P^ 

2 

Foe-^'t 


1 


P + a 

P + a 

3 

Fq(i - e”^"^) 

Foa 

a 

P(P + a) 

P(P + a) 

k 

Foe“®’*^(l - 

Fol) 

b 


■ (P +-a)(P + a + t) 

(P + a)(P + a + b) 

5 

Fq Bln at 

Foa 

a 


p2 + a2 

P2 + a2 

• 
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